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ABSTRACT 

An  experimental  study  of  airblast  atomization  was 
conducted  using  an  especially  designed  atomizer  in 
which  the  liquid  first  impinges  on  a  splash  plate,  then 
is  directed  radially  outward  and  is  atomized  by  the  air 
passing  through  two  concentric,  vaned  swirlers  that 
swirl  the  air  in  opposite  directions.  The  effect  oF 
flow  conditions,  air  mass  velocity  (mass  flow  rate  per 
unit  area,  paUa)  and  liquid  to  air  ratio  on  the  mean 
drop  size  was  studied.  Seven  different  ethanol  solu¬ 
tions  were  used  to  simulate  changes  in  fuel  physical 
properties.  The  range  of  atomizing  air  velocities  was 
from  30  to  80  m/s.  The  mean  drop  diameter  was  meas¬ 
ured  at  ambient  temperature  (295  K)  and  atmospheric 
pressure. 

NOMENCLATURE 

Dp  prefilmer  diameter,  m 

PE  peak  of  the  weight  distribution 
Re  Reynolds  number,  oDU/v 

SMD  Sauter  mean  diameter,  m 

U  velocity,  m/s 

W  mass  flow  rate,  kg/s 

W1  width  of  the  weiqht  distribution 
We  Weber  number,  DoU^/o 

p  density,  kg/m^ 

r  gamma  function 

a  surface  tension,  kg/s? 

u  dynamic  viscosity,  kg/ms 

a  efficiency  factor 


Subscripts-. 


INTRODUCTION 


The  fuel  spray  characteristics  have  a  great  influ¬ 
ence  on  the  performance  of  gas  turbine  combustors.  A 
change  in  the  flow  corlitions  or  in  the  physical  prop¬ 
erties  of  the  fuel  produces  a  change  in  the  spray  char¬ 
acteristics  of  the  fuel  injector.  The  latter  will  have 
a  gi  ater  influence  in  the  future  when  the  supply  of 
high-quality  fuels  cannot  be  satisfied  and  fuels  with 
different  physical  properties  must  be  used.  During  the 
past  few  oecades  many  reseachers  have  studied  the  ef¬ 
fect  of  liquid  properties  and  flow  conditions  on  atomi¬ 
zation  and  found  empirical  equations  for  the  type  of 
fuel  injector  investigated. 

The  work  of  Radcliffe  U)  showed  that  for  a  swirl 
atomizer  the  degree  of  atomization  depends  on  the  vis¬ 
cosity,  surface  tension,  mass  flow  rate  and  pressure 
drop  of  the  fuel.  He  found  the  exponents  0.6,  0.2, 
0.25,  and  -0.4  for  the  surface  tension,  viscosity, 
mass  flow  rate,  and  pressure  drop,  respectively. 

Jasuja  [2)  worked  on  fuels  having  surface  tension  rang¬ 
ing  in  viscosity  from  1.0xl0~6  to  93.0xl0~6  m2/s  find¬ 
ing  a  different  power  for  the  viscosity,  0.16.  He  ob¬ 
served  the  same  power  for  the  surface  tension,  0.6,  but 
the  variation  was  only  20  percent  and  was  accompanied 
by  a  large  variation  in  viscosity. 

Simmons  and  Harding  Q)  studied  the  atomizing  per¬ 
formance  of  six  simplex  pressure-atomizers  using  water 
and  kerosene,  liquids  with  almost  the  same  viscosity,  a 
30  percent  difference  in  density  and  a  water  surface 
tension  three  times  higher.  They  concluded  that  any 
difference  in  Sauter  Mean  Diameter  (SMD)  was  due  to  the 
difference  in  surface  tension  rather  than  density.  It 
was  found  that  the  power  for  the  surface  tension  is 
0.16  for  a  constant  liquid  pressure  and  0.19  if  the 
mass  flow  rate  was  held  constant. 

Merrington  and  Richardson  (4)  found  that  the  SMD 
for  a  plain-orifice  atomizer  was  proportional  to  the 
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viscosity  of  the  fuel  raised  to  the  0.2  power  and  in¬ 
versely  proportional  to  the  fuel  velocity,  for  liquid 
jets  injected  cross  stream  from  simple  orifices  into 
axial-flow  airstreams,  Ingebo  (5)  found  that  the  mean 
drop  diameter  was  proportional  To  the  product  of  the 
Weber  and  Reynolds  numbers  (WeRe)  raised  to  the  0.25 
power  for  WeRe  <  10“  and  proportional  to  the  0.4 
power  for  WeRe  >  106. 

Fraser,  Dombrowski  and  Routley  (6)  studied  the 
rotary  atomizer.  Their  studies  showetJ  that  the  SMD 
is  a  combination  of  a  constant  plus  a  term  including 
the  effects  of  surface  tension,  kinematic  viscosity, 
mass  flow  rates  and  relative  velocity  between  the  air 
and  the  fuel.  The  powers  0.5  and  0.21  were  found  for 
surface  tension  and  kinematic  viscosity,  respectively. 
Using  a  large  range  of  disc  types,  Friedman,  Gluckert 
and  Marshall  (2)  correlated  their  results  for  SMD  in 
terms  of  the  operating  and  liquid  variables  in  dimen¬ 
sionless  groups.  These  groups  present  the  viscosity 
raised  to  the  0.2  power  and  the  surface  tension  and 
density  raised  to  the  0.1  power. 

After  studying  the  experimental  data  on  prefilming 
types  of  airblast  atomizers  Lefebvre  (8)  concluded  that 
for  liquids  of  low  viscosity  the  main  Tactors  governing 
SMD  are  liquid  surface  tension,  air  density  and  air 
velocity,  whereas  for  liquids  of  high  viscosity,  the 
SMD  is  more  dependent  on  the  liquid  properties, 
especially  viscosity.  This  fact  had  been  observed  by 
Nukiyama  and  Tanasawa  (9) ,  Kim  and  Marshall  ( 10), 
Lorenzetto  and  Lefebvre  (11),  Rizkala  and  Lefebvre 
(2£) ,  El-Shanawany  and  LeTebvre  (_13),  and  Jasuja  (2). 
They  expressed  the  SMD  as  the  sum  of  two  terms,  the 
first  dominated  by  air  density  and  velocity,  and  the 
second  by  liquid  viscosity. 

An  experimental  investigation  was  conducted  to 
study  the  effect  of  mass  velocity  (mass  flow  rate/unit 
area,  o/vUa),  liquid  to  air  ratio,  and  liquid  properties 
on  the  spray  characteristics  of  two  fuel  injector 
modules  designed  for  high  temperature  and  high  pressure 
combustors.  The  experiment  was  conducted  in  an  open 
duct  facility.  The  SMD  of  the  spray  was  measured  at 
ambient  temperature  (295  K)  and  atmospheric  pressure. 
The  liquids  used  were  different  aqueous  solutions  of 
ethanol . 

The  measured  SMD  was  plotted  against  the  air 
mass  velocity,  liquid  to  air  ratio,  and  prefilmer  diam¬ 
eter  for  both  fuel  module  injectors.  The  data  were 
correlated  using  a  basic  equation  derived  by  Lefebvre 
(8). 

APPARATUS  AND  PROCEDURE 
Test  Facility 

A  schematic  of  the  test  facility  is  shown  in 
Fiq.  1,  and  a  photo  in  Fig.  2.  The  fuel  injector 
module  was  mounted  on  a  0.635-cm-thick  plate  and  ins¬ 
talled  on  the  end  of  a  15.25-cm-diameter  pipe.  Air  was 
supplied  by  the  Lewis  Research  Center's  air  system  with 
a  range  of  test  flow  rate  from  0.0169  tn  0,0684  kg/sec 
and  a  maximum  pressure  of  1.171x10“  Pa  at  the  fuel 
module  location.  A  pitot  tube  was  located  about 
22.5  cm  upstream  of  the  injector  and  connected  to  a 
manometer  board  which  was  used  to  set  the  pressure 
differential  across  the  module. 

A  pressurized  tank,  a  15.25-cm-diameter  schedule 
40  stainless  steel  pipe  1.22  m  long,  was  used  to  supply 
the  ethanol  solutions  to  the  fuel  injector  modules.  A 
pressure  regulator  was  used  to  keep  the  pressure  of  the 
nitrogen  in  the  tank  constant  at  5.15xl(F  Pa.  The 
liquid  mass  flow  rate  was  measured  usinq  a  rotameter 
previously  calibrated  for  each  solution  at  the  appro¬ 
priate  working  temperature. 


The  exhaust  and  liquid  collection  systems  con¬ 
sisted  of  a  0.16-cm-thick  stainless  steel  duct  about 
38.7  cm  in  diameter  and  3.8  m  long,  an  air  operated 
flow  amplifier  was  used  to  increase  the  velocity  of 
the  exhaust,  a  water  spray  was  installed  in  the  duct 
to  dilute  the  ethanol  solutions  and  avoid  any  flam¬ 
mable  mixture,  and  a  200  L  tank  was  used  to  collect 
the  solutions. 

Fuel  Injector  Modules 

The  fuel  injector  module  designs  used  in  these 
tests  are  shown  in  Fig.  3.  Each  fuel  injector  module 
consisted  of  two  concentric,  vaned  air-swirlers  that 
swirl  the  air  in  opposite  directions  to  create  a  zone 
of  high  shearing  action.  All  vanes  were  at  an  angle 
of  45  to  the  axial  direction.  The  liquid  was  supplied 
to  each  module  by  a  tube  located  in  the  central  cavity 
of  each  module.  The  liquid  flows  from  the  fuel  tube 
through  an  0.084-cm-diameter  discharge  opening  and 
impinges  on  a  splash  plate  mounted  on  the  downstream 
face  of  each  module.  This  splash  plate  breaks  up  the 
fuel  jet  and  directs  it  radially  outward,  where  the 
fuel  is  further  atomized  by  the  air  passing  through 
the  inner  swirler.  Additional  fuel  atomization  occurs 
in  the  shearing  region  between  flows  exiting  the 
counter  rotating  air  swirlers. 

Drop  Size  Measurements 

Drop  sizes  were  measured  using  a  Malvern  S.T. 

1800  Particle  and  Droplet  Size  Distribution  Analyzer. 
The  Malvern  instrument  is  a  nonintrusive  optical  system 
based  on  the  Fraunhofer  diffraction  of  a  parallel  mono¬ 
chromatic  light  beam  scattered  by  moving  droplets.  The 
transmitter  portion  of  the  Malvern  instrument  houses 
the  2-mW  helium-neon  laser  and  beam  expander,  which 
emits  an  approximately  9-mm-diameter  beam.  The  re¬ 
ceiver  consists  of  a  focusing  lens  (Fourier  transform 
lens),  a  multielement  photoelectric  detector,  beam 
alinement  knobs,  lamps,  and  an  indicator.  A  computer 
with  an  8  K  memory  receives,  stores,  and  processes  data 
inputs  from  the  detector.  A  teletype  with  a  hard  copy 
printer  is  used  for  data  output.  The  output  is  dis¬ 
cussed  in  the  appendix.  Two  data  points  were  taken  at 
each  condition  and  stored  in  the  computer  memory. 
Measurements  were  made  at  the  center  line  of  the  spray 
at  a  distance  of  7.62  cm  downstream  of  the  fuel  injec¬ 
tor  module. 

RESULTS  AND  DISCUSSION 
Mass  Velocity 

The  effect  of  mass  velocity  (mass  flow  rate  per 
unit  area,  oaua)  on  for  the  fuel  inJector 

modules  investigated  is  clearly  shown  in  Figs.  4  and  5. 
These  figures  show  in  general  that  the  SMD  decreases 
with  an  increase  in  the  mass  velocity.  The  same  effect 
was  observed  for  both  modules.  Changes  in  mass  veloc¬ 
ity  were  obtained  by  changing  the  total  air  flow  rate 
through  the  injector  modules  while  keeping  the  avail¬ 
able  flow  area  (including  both  swirlers)  constant. 

The  range  of  mass  velocities  was  from  37  to  112  kg/ms 
and  the  calculated  flow  area  was  6.428xl0-4  m?  for 
injector  module  1  and  4.519xl0~4  m?  for  injector 
module  2.  The  area  was  calculated  at  the  upstream  side 
of  the  injector  modules. 

Liquid  To  Air  Ratio 

Tests  were  conducted  to  study  the  effect  of  liquid 
to  air  ratio  on  SMD.  These  tests  covered  a  range  of 
liquid  to  air  ratios  from  0.0147  to  0.0462  for  module  1 
and  from  0.0202  to  0.0636  for  module  2.  Figures  6  and 
7  show  no  effect  of  the  liquid  flow  rate  on  SMD  when 
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the  air  velocity  rate  is  kept  constant  at  different 
concentrations  of  ethanol. 

Liquid  Properties 

Many  tests  were  made  to  determine  the  atomizing 
performance  of  two  especially  designed  fuel  injector 
modules  operating  with  liquids  of  different  physical 
properties.  A  number  of  aqueous  solutions  of  ethanol 
were  prepared  representing  the  following  range  of 
liquid  properties: 

surface  tension  =  0.0290  to  0.0555  kg/s^ 

dynamic  viscosity  =  0.001226  to  0.002684  kg/ms 

density  »  890.7  to  988.0  kq/m^ 

Samples  of  the  solutions  were  analyzed  using  standard 
laboratory  techniques  to  measure  surface  tension,  vis¬ 
cosity  and  density.  Table  !  presents  the  results  of 
those  measurements. 

The  effect  of  ethanol  concentration  on  SMD  for 
the  fuel  modules  studied  is  shown  in  Fig.  8.  Both 
figures  show  a  decrease  in  the  SM0  with  an  increase 
in  the  ethanol  concentration,  i.e.,  decrease  in  sur¬ 
face  tension. 

Figure  8  shows  a  comparison  between  data  from 
Ref.  12  for  an  airblast  atomizer  and  data  from  this 
investigation.  Both  fuel  injector  modules  produced 
smaller  droplets  than  the  airblast  atomizer  for  the 
same  air  velocity,  60  m/s.  The  only  difference  is  the 
liquid  flow  rate,  but  Figs.  6  and  7  showed  the  SMD  is 
not  affected  by  changes  in  the  liquid  to  air  ratio. 

This  figure  shows  the  benefits  created  by  the  high 
shearing  action  between  the  flows  exiting  the  counter 
rotating  air  swirlers  compared  with  the  airblast  atom- 
izer  of  Ref.  12.  Note  that  the  liquid  surface  tension 
was  used  as  a  parameter  because  the  liquid  used  in 
Ref.  12,  aqueous  solutions  of  Butan-2-ol,  have  almost 
the  same  physical  properties. 

Linear  Scale 

Two  identically  designed  fuel  injector  modules 
were  used  in  this  investigation.  The  only  difference 
was  the  size  and  the  number  of  swirler  vanes  (Fig.  3). 
Module  2  is  approximately  20  percent  smaller  in  diam¬ 
eter  than  module  1,  having  a  prefilmer  diameter.  Dp, 
of  1.2  and  1.5  cm,  respectively.  The  influence  of 
atomizer  scale  on  SMD  is  illustrated  in  Fig.  9.  This 
figure  shows  module  2  producing  smaller  droplets  than 
module  1  under  the  same  operating  conditions. 

Data  Analysis 

The  experimental  data  gathered  in  this  investiga¬ 
tion  were  correlated  usinq  the  basic  equation  derived 
by  Lefebvre  (8)  with  the  experimental  constants  of 
Ref.  13 


As  stated  in  Ref.  13,  the  experimental  constants  0.073 
and  0.015  may  have  to  be  modified  by  an  efficiency 
factor,  whose  values  will  depend  on  the  atomizer 
design  and  will  take  into  account  the  presence  of  ex¬ 
traneous  devices,  such  as  air  swirlers,  and  the  differ¬ 
ent  methods  of  drop  size  measurement.  Figures  10  and 
11  compare  the  SMD  measured  in  the  present  investiga¬ 
tion  with  values  predicted  by  Eq.  (1).  Good  agreement 
is  shown  in  Fig.  11  for  a  value  of  9  of  0.59. 


Figure  10  shows  good  agreement  in  the  high  velocity 
region  for  the  same  value  of  9,  but  does  not  describe 
the  experimental  data  very  well  in  the  intermediate  to 
low  velocity  region  i.e.,  less  than  60  m/s,  correspond¬ 
ing  to  drop  sizes  greater  than  60  u. 

SUMMARY  AND  CONCLUSIONS 


An  experiment  was  conducted  at  atmospheric  pres¬ 
sure  to  determine  the  effect  of  liquid  physical  prop¬ 
erties  and  flow  conditions  on  the  Sauter  Mean  Diameter, 
SMD,  using  two  geometrical ly-simi lar  research  airblast 
atomizers  designed  for  high  pressure  and  temperature 
combustors. 

After  studying  the  effects  of  the  different  vari¬ 
ables  involved  in  this  investigation  it  is  found  that: 

1.  The  SMD  of  the  spray  decreases  with  increases 
in  ethanol  concentration  due  to  changes  in  the  physical 
properties  as  shown  in  Table  I. 

2.  Increasing  the  air  velocity  decreases  the  SMD 
which  varied  inversely  with  air  velocity,  of  all  the 
variables,  air  velocity  has  the  most  dominant  effect  on 
the  atomization  process  of  the  fuel  injector  modules 
investigated. 

3.  An  increase  in  atomizer  scale  increases  the 
mean  drop  diameter. 

4.  The  air  to  liquid  ratio  has  no  measurable  ef¬ 
fect  on  the  SMD  of  sprays  produced  by  any  of  the  two 
fuel  injector  modules. 

5.  The  SMD  performance  of  the  airblast  atom¬ 
izers,  when  spraying  in  stagnant  air  at  atmospheric 
pressure,  is  predicted  with  a  reasonable  degree  of 
accuracy  by  the  correlation: 


SMD  = 


0.1237, 


0.5 


0.02541 


(•5) 


(2) 


for  the  following  ranqe  of  test  conditions 
Surface  tension  =  0.0290  to  0.0555  kq/s? 

Dynamic  viscosity  =  0.001226  to  0.002684  kg/ms 
Liquid  density  =  890.7  to  988.0  kg/m^ 

Air  velocity  =  30  to  BO  m/s 

Liquid  to  air  ratio  =  0.0147  to  0.0776 

Air  density  was  not  varied  appreciably. 
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APPENDIX  -  OUTPUT  FROM  THE  MALVERN  S.T.  1800 
PARTICLE  AND  DROPLET  SIZE  DISTRIBUTION  ANALYZER 

The  Malvern  instrument  is  a  nonintrusive  optical 
system  based  on  the  laser  diffraction  principle.  This 
instrument  uses  the  Rosin-Rammler  weight  distribution 
model.  The  Rosin-Rammler  distribution  is  defined  as 
follows 

p(x)  -  exp  |-(x/PE)W'| 

where  P(x)  is  the  weight  or  volume  fraction  of  par¬ 
ticles  in  the  range  x  to  x  ♦  dx  where  x  is  in 
microns.  The  parameters  PE  and  W'  characterize  the 
peak  of  the  weight  distribution  and  its  width.  PE  is 
in  microns  and  W'  is  a  dimensionless  number  usually 
in  the  range  from  one  for  very  wide  weight  distribu¬ 
tions  to  10  for  very  narrow  weight  distributions.  The 
values  of  the  two  parameters,  PE  and  W1,  which  give 
the  minimum  error,  E,  define  the  size  distribution. 
Using  the  values  of  PE,  W'  and  the  gamma  function  the 
Sauter  Mean  Diameter  (SMD)  can  be  calculated.  The  fol¬ 
lowing  formula  is  used 


where  SMD  is  in  microns  and  r  is  the  tabulated 
gamma  function. 

Figure  13  shows  an  example  of  the  output.  The 
first  line  of  output  is  the  peak  PE,  width  W‘,  and 
error  E  of  the  distribution.  The  first  column  gives 
the  droplet  size  ranges  in  microns.  The  next  three 
columns  are  the  spray  distributions  as  percent  weight 
fraction,  cumulative  percent  by  weight,  and  normalized 
percent  by  number  density.  The  last  two  columns  are 
the  calculated  and  actually  measured  energy 
distributions. 
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TABLE  I.  -  LIQUID  PROPERTIES 


Level 

tested 

Solution, 

percent 

ethanol 

u» 

kg/ms 

°*  o 

kg/s2 

p’ 

kg/m-3 

1 

5 

0.001226 

0.0555 

988.0 

2 

7.5 

.001264 

.0530 

985.3 

3 

10 

.001438 

.0485 

980.8 

4 

20 

.001973 

.0395 

969.4 

5 

30 

.002488 

.0340 

952.4 

6 

40 

.002684 

.0325 

936.0 

7 

60 

.002424 

.0290 

890.7 

schematic.  (Not  to  scale. ) 
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Figure  10.  -  Comparison  of  measured  and  predicted  values  of 
SMD  from  equation  (2),  Module  1. 
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Figure  11.  -  Comparison  of  measured  and  predicted  values 
of  SMD  from  equation  (2),  Module  2. 


>  PE  *  +102  0  W  =  +2  4  E  =  00313752 


D  =  +562. 86  > 
D  -  +261. 71  > 
D  ■  +160. 29  > 
D  ■  +112. 86  > 
D  =  +84. 29  > 
D  ■  +64. 57  > 
D  -  +50. 29  > 
D  *  +38.86  > 
D  ■  +30. 29  > 
D  =  +23.71  > 
D-  +1°,  57  > 
D*  +14.57  > 
D  -  +11. 43  > 
D-  +9.14  > 
D  -  +7. 14  > 


+261. 71 
+160. 29 
+112.86 
+84. 29 
+64.57 
+50.29 
+38.86 
+30.29 
+23.71 
+18.57 
+14.57 
+11. 43 
+9.14 
+7.14 
+5.71 


+0.01% 
+5. 18% 
+22. 76% 
+25. 17% 
+18.50% 
+11.64% 
+7.34% 
+411% 
+2  31% 
+1. 31% 
+0.73% 
+0.41% 
+0. 22% 
+0.14% 
+0. 07% 


R  ■  +99. 99% 
R  •  +9481% 
R  ■  +72  05% 
R-  +46.88% 
R  -  +28. 38% 
R  -  +16.74% 
R  *  +9.39% 
R-  +5.28% 
R«  +2.97% 
R-  +1.66% 
R  «  +0.93% 
R-  +0.52% 
R  ■  +0.31% 
R-  +0.17% 
R-  +0.10% 


N-  +0.00% 
N-  +0.03% 
N  ■  +0.54% 
N-  +1.60% 
N  -  +2  73% 
N-  +3.74% 
N-  +5.05% 
N-  +6.06% 
N«  +7.14% 
N-  +8.42% 
N  *  +9.77% 
N-  +11.39% 
N  =  +1208% 
N-  +15.39% 
N-  +16.06% 


0  0698 
O  1005 
O  1380 
C  *  1707 
C -  1970 
O  2044 
C-  1980 
C  -  1783 
C  -  1542 
0  1303 
C  -  1067 
O  0867 
C  -  0691 
O  0541 
C  *  0419 


A  -  0756 
A  ■  1157 
A  -  1424 
A  ■  1646 
A  -  1780 
A  -  2002 
A  -  2047 
A  -  2002 
A  -  1869 
A  -  1557 
A  -  1112 
A  -  0623 
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A ■  0222 
A  -  0178 


Figure  12  -  Example  of  output  from  Malvern  particle  and  droplet  size  distribution  analizer. 
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